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ELECTRON  VELOCITY  SHEAR  INSTABILITY  IN  THE  AURORAL  IONOSPHERE 


I.  Introduction 

Electron  streams  flowing  in  magnetic  fields  are  a  common  occurrence  in 
terrestrial  and  astrophysical  plasmas,  e. g. ,  the  earth's  ionosphere  and 
magnetosphere  [Arnoldy ,  1974-  Anderson  and  Vondrak,  1975],  and  solar  flares 

[S  turrock,  1974].  These  streams  have  been  observed  flowing  parallel  or 
perpendicular  to  the  magnetic  field  with  the  flow  often  spatially  varying 
perpendicular  to  the  magnetic  field.  It  is  well  know  that  this  flow 
configuration  is  unstable  (Chandrasekhar ,  1957]  and  subject  to  a  Kelvln- 
Helnho It z-like  velocity  shear  instability.  For  the  ionosphere-magnetospheric 
plasma,  Webster  [1957]  suggested  the  possible  connection  between  small  scale 
structures  in  auroral  arcs  and  instabilities  arising  from  electron  x  B_  flows 
perpendicular^  to  the  magnetic  field.  However,  there  is  considerable  evidence 
[Arnoldy  1981,  and  references  therein]  that  strong  inhomogeneous  electron 
flows  parallel  to  the  geomagnetic  field  also  exist  near  auroral  arcs.  To  our 
knowledge,  the  instabilities  of  transverse  sheared  electron  streams  flowing 
parallel  to  the  magnetic  field  has  not  been  fully  discussed  in  a  space  plasma 
physics  context. 

In  the  following  we  present  a  linear  kinetic  theory  of  transverse 
velocity  sheared  electron  flows  parallel  to  the  magnetic  field  (Harrison, 
1963;  Harrison  and  Stringer_,  1963;  Mikhailovskii  and  Rukhadze ,  1966  Rome  and 
Briggs,  1972).  Our  treatment,  which  includes  density  gradients,  arbitrary 
Te/T^,  background  ions  and  electron  collisions,  is  applicable  to  the 
low  B  (particle  pressure/magnetic  pressure)  weakly  ionized  auroral  ionospheric 
plasma.  We  show  that  the  interface  between  downward  flowing  hot  electrons  and 
upward  flowing  return  current  cold  electrons  is  unstable  to  a  high  frequency 
electron  velocity  shear  driven  instability.  In  Section  II  we  derive  a  general 
Manuscript  submitted  April  13,  1982. 
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kinetic  theory  and  present  our  principal  results.  In  Section  III  we  summarize 
our  findings  and  discuss  the  nonlinear  implications  of  this  instability. 


II.  Theory 

A.  Physical  Mechanism  of  the  Instability 

We  first  present  a  simple  physical  picture  of  the  instability  driven 
by  a  transversely  sheared  electron  flow  along  a  magnetic  field  (Rome  and 
Briggs,  1972).  The  magnetic  field  B  =  B  e^  is  taken  to  be  constant  in  space 
and  time.  The  ions  are  assumed  to  form  an  immobile,  homogeneous  background, 
while  the  electrons  are  assumed  to  have  an  inhomogeneous  flow 

A 

velocity  V  =  V  (x)e  such  that  V  (x)  -  (x-x  )  3V  /3x  with  3V  /3x  <  0. 
oe  oe  z  oe  o  oe  oe 

We  impose  upon  this  configuration  a  two-dimensional,  perturbed  electric  field 

55  *  5Evew  +  <5E,e,  as  shown  in  Fig.  (1).  The  perturbation  force  in  the  z- 
y  y  z  z 

direction  acting  on  an  electron  fluid  element  at  x  **  xQ,  can  be  written 


<5F 


ez 


<5Fa  +  <5F 
ez 


d 

ez 


where 


5F 


ez 


-  e  6E 


<5F 


ez 


-  m  <5V  (3V  /3x) 

e  ex  oe 


and  <5V  3  "  c6E  /B  is  the  perturbed  6E  x  B  drift  due  to  5E  .  Here,  5F3  is 

ex  y  ~  ~  y  ez 


the 

"acceleration"  force 

due 

to  5E  which  acts  to  neutralize 
z 

the  charge 

imbalance,  while  6Fd  is 

ez 

the 

"deceleration"  force  due  to  the 

convection 

term 

(<SV  •  V)V  and  opposes 

fFa 

ez 

when  3V  /3x  <  0.  Thus,  when  $Fd 
oe  ez 

>  5Fa  ,  the 
ez* 

force  imbalance  on  an  electron  fluid  element  leads  to  an  enhancement  of  charge 
separation,  and  hence,  instability.  For  electrostatic  fluctuations. 
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<SEz  =  -  i  kz<5$  and  »  -  i  k  6<f>,  where  6<p  is  the  electrostatic  potential; 

this  gives  the  following  instability  criterion  |  ( l/£2  )3V  /3x|  >  |k  /k  |.  We 

e  eo  z  y 

prove  this  criterion  rigorously  in  the  next  sections. 


B.  Dispersion  Equation 

We  consider  a  plasma  configuration  as  shown  in  Figure  2.  Note  that 
the  ions  are  assumed  to  have  no  flow  in  the  z-direction  and  that  the 
temperature  is  taken  to  be  constant.  We  consider  only  electrostatic 
oscillations  and  assume  perturbed  quantities  vary  as  exp[i  (k  y+k  z-ut)] 

y  2 

with  k2/k2  «  l.  We  assume  a  weakly  inhomogeneous  plasma  (r2^  «  L2  where 
r  is  the  mean  ion  Laraor  radius  and  L  is  the  scale  length  of  the  density  and 

1-iX 

velocity  gradients)  and  make  use  of  the  local  approximation 

(k2  >>  32/ 3x2  »  i/L2).  We  assume  that  ft?  <<  u>2  <<  ft2  so  that  the  ions  can  be 
y  l  e 

considered  unmagnetized  and  the  electrons  magnetized.  Finally,  we  include 
electron  collisions  in  the  analysis,  i.e. ,  electron-ion  collisions. 

The  unperturbed  electron  distribution  function  is  given  by 


Foe(vi’  V  X)  "  noe(X)(lTVer3/2  e*P[-<v^Virvoe0O)2)/v2] 


where  v2  =  v2  -v2  v„  *  v  and  X  =  x  -  v  /ft  are  constants  of  motion. 
1  x  y  a  z  y  e 

1/2 

Here,  v  =  (2T  /m  )  and  ft  *  eB  /m  c  are  the  electron  thermal  velocity  and 
eee  eoe 

cyclotron  frequency,  respectively.  We  expand  (1)  about  x  =  0  and  obtain 


v  3  F 

F  (v2  v  X)  -  Fm  (v2  v  )  -  — y -2®. 

oe^  i’  r  oe  '  1’  r  ft  3x 
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where 


Fm  (v2  v  )  =  n  (itv2)  exp [ - ( v 2+( v  -V  ) 2 ) / v2  I 

oe  1  I  oe  e  1  II  oe  e 


x  =  o 


(3) 


and 


3  F  I 

— °£  „  Fn 

3x  1  x=0  oe 


0 

The  equilibrium  electron  drifts  associated  with  (2)  are  Vqx  =  0, 

Ve  »  -  (v2/3Q  )3  in  n  /3x  and  V6  =  V  where  Ve  **  f  d3v  v  F 
oy  e  e  oe  oz  oe  J  ~  oe 

The  electron  equilibrium  pressure  tensor  is  defined  as 

P  **  f  d3v  (v  -  V  8)(v  -  V  6)  F  and  is  given  by 
*oe  J  ~  ~  ■'Q  oe  °  7 


3  in  n 


oe 


3x 


2(v  -V  )V  3  in  V 
_  J _ oe_  oe  _  oe_ 

?  3x 


(4) 


x=0 


P 


oe 


=n  T 


oe  e 


(5) 


where  A  *  -  (l/f2  )3V  /3x.  It  is  the  anisotropy  in  the  pressure  tensor  that 

e  e  oe 

provides  the  free  energy  to  drive  the  instability. 

The  unperturbed  ion  distribution  function  is  assumed  to  be 


"  noi(^Vi) 


-3/2 


exp 


[-(V* 


+(Vvdi>2 


+v‘)/v2) 

Z  X 


(6) 


where  v^  =  (2Tj/m^)^2  is  the  ion  thermal  velocity,  =  (v2/2Sh)  3  in  nQ^/3x 
is  the  ion  diamagnetic  drift  velocity  and  =  eB^/n^c  is  the  ion  clclotron 
frequency. 

The  electrostatic  dispersion  equation  is  obtained  from  Poisson's  equation 


V25<S  =  47r(6n  -  Sn  ) 

e  i 


(7) 
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where  $<J>  is  the  perturbed  potential,  5n^  is  the  perturbed  density  of 
species  o,  i.e,  *  /  d3v  df^,  arid  <5fQ  is  the  perturbed  distribution 
function  of  species  o.  We  define  *  -  (Aire^/k2  )6nQ  so  that  the  dispersion 
equation  is 


D(oj,k)  =  1  +  +  Xe  =  0 


(8) 


We  obtain  xe  by  making  use  of  the  BGK  or  Krook  model  to  describe  the 
effects  of  electron  collisions  (Kadomtsev,  1965).  That  is,  we  consider 

(v  +  v  •  - -e-  --  v  x  B  •  )<Sf  =  -  —  ik*  <5<t>-v  (<Sf  -  -~e  F  ).  (9) 

v3t  -  3x  m  c  -  ~  3vJ  e  n  -  3v  e  e  n  oe 

-  e  ~  e  o 


where  represents  the  electron  collision  frequency.  Following  standard 

techniques,  i.e.  ,  the  method  of  characteristics,  we  determine  Sf  from  Eqs. 

e 

(2)  and  (9);  then  integrating  over  velocity  space  we  finally  obtain 


X,  = 


2oo2  k  V ,  .  k  , 

_.P±n-r’-  -  2  -d±)7(r  )  -  1  _£  1 

,2  2  X '  ’  ^e  k  v  2  k  P. 

k^v4  z  e 


dV 


z  e 


oe 

dx 


z  (c  )]r  Cb  )]  X 
e  o  e 


IV 

[i+  t— z  (c  )  r  (b  )  ] 

kv  e  o  e  1 

z  e 


-1 


(10) 


where  C  =  (u  +  iv  -k  V  )/k  v  ,  u2  =  4irne2/m  ,  T  (b  )  =  exp(-b  )  I  (b  ), 
e  e  z  oe  ze  pe  e  o  e  e  o  e 

b  =  k2r2  /2,  rT  =  v  /fl  ,  I  is  the  modified  Bessel  function  of  order 
e  y  Le  Le  e  e  o 

0,  V,  =  -  (v2  /2S2  )  3  £n  n/3x,  Z  is  the  plasma  dispersion  function 
de  e  e 

and  Z'(c)  =  dZ/d;.  The  ion  response  is  simply  given  by 


2ui2. 

— ^  U  +  C.  Z  (c.)J 

k2v2  1  1 


(11) 
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where  <o2  =  4Trne2/m. ,  C,  «  (w  -  k  V  )/kv,  and  V,.  *=  (v2/2P,  )3  Jtn  n/3x.  Thus, 

pi  11  y  ai  i  ai  1  i 

Eqs.  (8),  (10)  and  (11)  describe  electrostatic  waves  In  the  presence  of  a 

sheared  parallel  electron  velocity  for  the  frequency  regime  P2  «  id2  «  P2. 


C.  Analytical  Analysis 

In  order  to  gain  insight  into  the  nature  of  the  electron  velocity  shear 

instability,  we  analytically  solve  Eq.  (8)  in  two  limiting  cases:  "cold" 

electrons  (Tg  <<  Tj )  and  "hot"  electron  (Te  >>  T^).  For  simplicity,  we  expand 

about  x  =  xQ  and  an  consider  electron  velocity  profile  such  that 

V  =  (x-x  )  (3V  /3x)  so  that  V  (x=x  )  =  0. 

oe  o  oe  x=x  oe  o 

o 

1.  Cold  electrons  (Tg  <<  ) 

In  this  limit,  the  contribution  from  the  ions  can  be  ignored 
since  <<  1  and  the  dispersion  equation  can  be  written 


D(m,k)  =  1  +  x  =0 


(12) 


In  order  to  simplify  Eq.  (10)  we  make  the  following  assumptions: 

(to  +  i  v  )/k  v  <<  1  and  k2r2  <<  1.  For  these  conditions  we  note  that 
e  z  e  y  Le 

Z(c  )  =  -  1/s  -  l/2s2,  Z'U  )  =  1/C2  and  T  (b  )  =  i  -  k2r2  12.  We  find  that 

e  e  e  e  e  o  e  y  Le 


..,2 


xe=l 


_  PS. 

P2 

e 


pe  z  r  z  1 
k  ^k  +  P 


(u+iv  )2  y  y 
e 


dV 

i 

Tx 


0£j 


w  ai  i  v  , 

_  Pe_  P.£  _  °  In  n.  r  __  e_  \-l 

to+iv  P  k  3x  *•  u+iv  ' 
c  e  e 


(13) 


Making  use  of  Eqs.  (12)  and  (13)  we  arrive  at  the  following  equation 
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(14) 


D(w,k)  *  (l+o2)oi2  -  (ivg-  cr/k^Ln)oj  -  A  ■  0 


where  to  =  (tiH-iv  )/to  ,  v  ■  v  /m  ,  a  =  u  /Q  ,  L  =  (3  £n  n/3x)  * 
e  pe*  e  e  pe  pe  e  n  * 

/  / 

A  ■  (k  /k  )(k  /k  +  V  /SI  ),  and  V  =  3V  /3x.  The  solution  of  Eq.  (14)  is 
z  y  z  y  oe  e  oe  oe  n 


w  = -|  [(ive~a/kyLn)^{(ive-a/kyLn)2  +  4A(  1+a2)} ^2 ] ( l+o2)  1 


(15) 


We  further  simplify  Eq.  (15)  by  examining  the  following  limits. 


a.  Collisionless  plasma  (v  =  0) 


In  the  collisionless  limit  (v  =  0),  Eq.  (15)  reduces  to 

e 


[-  t-V  ±  {-  ~  +  4A( l+o2)}  1/2 ]  (i  +  a2)-1 

kyLn  k2L2 
J  y  n 


Instability  can  occur  when 


(16) 


A 


3V  . 

_°* )  <  -  I 

3v  >  S  L 


1+a2  k2L2 
y  n 


(17) 


In  the  limit  L  ■*  «,  this  instability  criterion  is  identical  with  that  derived 
n 

from  the  simple  physical  picture  presented  in  Sec.  IIA.  The  important  aspect 
of  Eq.  (17)  is  that  the  density  gradient  is  a  stabilizing  influence 
(Hikhailovskii  and  Rukhadze,  1966).  Moreover,  the  density  gradient  is  most 
effective  in  stabilizing  the  instability  when  or2  »  1  (i.e.,  w2e  »  fi2)  and 
kyLn  <<  1.  The  latter  condition  indicates  that  the  density  gradient  acts  to 
stabilize  long  wavelength  modes  before  short  wavelength  modes,  so  that  a  "long 
wavelength  cut  off"  should  exist  in  inhomogeneous  plasmas.  We  also  note  that 
the  instability  attains  maximum  growth  for 
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(18) 


k  ,  3V 

_z  „  1  _oe 

k  £2  3x 

y  e 


b.  Collisional  plasma  (v  *  0) 

e 


We  consider  a  collisional  plasma  (vg  *  0),  but  neglect  the  density 
gradient  for  simplicity  (L^  +  0).  In  this  limit,  Eq.  (15)  becomes 


S  =  i  /  (1  +  a2)'1  [1  ±  (1  -  4-A(l^i)_)l/2] 


(19) 


A  necessary  condition  for  instability  is 


3V 


oe 


3x 


)  <  0 


(20) 


which  is  the  same  condition  as  Eq 
by 


(17)  for  L 

n 


The  growth  rate  is  given 


«  r_L'l±  i 
;  ->J 
l+a2 


1/2 


^e  1+2  a  2 
2  l+a2 


for  v2  <<  4A(1  +  a2) 
e 


(21) 


and 

y  =  JAL  _  _  ®L  v  for  v2  »  4A(  1  +  a2)  (22) 

.  .  j  e  e 

v  1+a^ 
e 

Equation  (21)  is  the  weakly  collisional  limit  and  Eq.  (22)  is  the  strongly 
collisional  limit.  A  key  parameter  in  determining  the  strength  of  the 
instability  in  a  collisional  plasma  is  a  =  “'pg/^g*  In  fact,  for  a2  »  1  it 
is  found  that  collisions  prevent  the  mode  from  becoming  unstable.  However, 
for  a2  -*■  0,  instability  presists,  even  for  a  large  collision  frequency. 
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2.  Hot  electrons  (Tfi  »  T^) 

In  this  limit  the  contribution  from  the  ions  cannot  be  ignored  and 
the  dispersion  equation  is 


D(w,k)  -  1  +  x±  +  Xe 


(23) 


Again,  to  simplify  Eq.  (23)  and  to  focus  on  the  role  of  the  ions,  we  make  the 
following  simplifying  assumptions:  =0,  -*•  00 ,  m  »  ^2ve»  ^yrLe  «  * 

and  u  »  k-zv^.  The  susceptibilities  xq  are  now  given  by 

oj 2  k 
_pe  _z 

2  k 

u*  y 

and 


k  ,  3V 

r  z  +  i-  _  22.) 

D.  Sv  j 


(24) 


e 


xi 


so  that  the  dispersion  equation  is 


m2  u2 

l  +  .£1  _  _P£ 

ft2  OJ2 


n 


3V 


oe 

3x 


m 


or 


m 

(1  +  a2)  cu 2  —  (A  +  — 6 )  =  0 

m. 

l 


Clearly,  instability  can  occur  when 


A 


3V 


<  - 


n 

e 


(25) 


(26) 


(27) 


(28) 
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Increased  stability  is  found  when  finite  electron  temperature  is  taken  into 


account  [Harrison,  1963 1 .  The  influence  of  the  ions  is  to  make  the 

instability  criterion  more  stringent  than  the  cold  electron  case, 
i.e.  ,  A  <  0.  Physically,  this  is  due  to  the  fact  that  the  ions  are  no  longer 
immobile  and  are  able  to  respond  to  the  perturbation  field  6Ez  (see  Fig.  1) 
Also,  Eq.  (28)  indicates  that  the  instability  is  easier  to  excite  in  "heavy 
ion"  plasmas. 

However,  we  can  show  that  an  instability  can  be  excited  even  when  Eq. 

(28)  is  not  satisfied,  i.e  ,  A  >  -  m  /m. .  In  this  case  it  is  crucial  that 

e  i 

electron  resonance  terms  be  retained  in  Eq.  (24),  which  have  been  neglected. 

\  1 0 

We  expand  the  electron  Z  function  as  Z (?  )  =  -  1/c  -  1/t2  +  iff1'  exp(-r2) 

e  e  e  e 

*  1  j 

and  Z  (4e)  =  1/ 42  “  ir  2?g  exp  (-C2)  so  that  dispersion  equation  becomes 


m  .  2wz 

D(w,k)  *  1  +  a2 - (A  +  —  )  +  iff*  ^  c  exp(-£2)  A  « 

in.  ,996  £ 


•9 

u>- 


X2v2 

z  e 


(30) 


where  £  =  w/k  v  .  We  assume  A  +  m  /m,  >  0  so  that 

e  z  e  e  i 


m 

r,  .  1  .,1/2 

w  =  w.,  [1  +  —  AJ 

r  £h  m 


(31) 


1/2 

where  =  u  j,/(l  +  o2)  .  Assuming  y  <<  it  is  easily  shown  that 


Y  =  -  ff*^  — -P—  X,  exp(~c2)  A  (1+a2)  *  to 

k2v2  6  6 

z  e 


(32) 


so  that  instablity  can  occur  when  A  <  0;  the  same  instability  criterion  as  in 
the  cold  electron  case.  An  important  aspect  of  this  instability  is  that  it 
has  a  real  frequency  associated  with  it  (Eq.  (31)),  while  the  other 
instabilities  discussed  are  purely  growing  (i.e.,  u>  =  0)  in  the  absence  of  a 


density  gradient.  We  add  that,  as  in  the  cold  electron  limit,  both  a  density 
gradient  and  electron  collisions  have  a  stabilizing  influence  on  this  resonant 
electron  velocity  shear  instability. 


D.  Numerical  Results 


To  better  illus crate  the  various  limits  presented  in  the  previous 

subsection,  we  present  a  set  of  curves  which  solve  Eq.  (8)  exactly  for  a  broad 

range  of  parameters.  The  important  parameters  that  are  varied  are 

T  /Tj  ,  k  rT  ,  r.  /L  ,  V  /ft  L,  (3 V  /3x  =  V  /L),  v  /Q  and  a  /ft  .  -  We  also 
e'i*  yLe  Le  n  oe  oe  o  e  e  pe  e 

consider  an  0+  plasma  for  these  cases,  i.e.,  m^  =  16  mp.  For  applications  to 
the  auroral  ionosphere  we  use  parameters  typical  of  regions  in  and  near 
discrete  auroral  arcs.  The  transverse  velocity  shear  for  electron  flow 


parallel  to  the  geomagnetic  field  can  be  approximated 

by  | 3V  /3x |  «  |V  /Li  where  L  is  a  typical  latitudinal  arc  width 

oe  o 

and  V  =  vL  -  v  is  c.-.e  difference  in  parallel  streaming  velocities  between 
o  h  c 

the  hot  downward  flcving  current  carrying  electrons  (vh)  and  the  upward 
flowing  cold  retcrr.  current  electrons  (vc).  Typical  values  for  the 
latitudinal  width  of  discrete  arcs  are  in  the  range  L  =  100  m  —  10  km  [Davis , 
19781  .  For  we  take  i„  2  en  v  -  10  A/m  [Anderson  and  Vondrak,  1975 

*  C  il  ,  C  C  C  - - 

A  5 

and  references  therein]  with  nc  =  10  -  10  cm  giving  vc  =  102  -  103 

cm/sec.  For  v.^  we  consider  the  field-aligned  downward  flowing  electrons  with 

9 

energies  E  =  1  -  10  keV  [ Arnoldy ,  1981]  giving  =  10  cra/sec.  As  a  result, 

4  -1 

an  estimate  for  the  transverse  shear  is  1 3V  /3x|  =  |  (v  -  v  )/L|  =  10  sec 

oe  he 

—2  6  — | 

with  L  =  1  km  giving  for  the  parameter  iV  /ft  L|  =  10  where  ft  =10  sec 

0  0  0 

In  Fig.  3  we  plot  w/u  vs.  T  /T..  The  solid  lines  show  the  growth 

pe  e  l 

rate  y/w  while  the  dashed  lines  show  the  real  frequency  u)  /u)  .  The 

pe  r  pe 

parameters  used  are  k  r.  *  0.1,  m  l  Si  =  0.5,  v  /w  *  0, 
r  y  Le  pe  e  e  pe 


Fig.  3  —  Plot  of  cj/cjpg  vs.  Tg/Tj  for  strong  shear  (V0/f2eL  =  —  0.02)  and  weaker 
shear  (V0/J2eL  =  —  0.01).  Note  stabilization  for  Te /T;  —  1  for  weak  shear  case. 


/ 


r,  /L  -  1.0  x  10  ^  for  V  /0  L  ■  -  0.01  we  choose  k_/k„  -  0.005,  and 
Le  n  o  e  z  y 

for  V  /!)  L  *  -  0.02  we  choose  k_/k„  -  0.010. 
o  e  z  y 

We  first  discuss  the  strong  shear  case,  Vq/9.^L  **  -  0.02.  We  note  that 

the  instability  criteria  A  =  (k  /k  )  (k  /k  +  V  /!!  L)  <  0  and 

zyzy  oe 

A  =*.  (k  /k  )  (k  /k  +  V  /52  L)  <  -  m  /m.  are  satisfied  for  the  parameters 

zy  zyoe  ei 

chosen.  As  expected,  a  nonresonant  instability  occurs  for  the  entire  range  of 
Tg/T^  considered.  Note  that  the  growth  rate  is  slightly  smaller  for  Te/T^  ■ 
10^  than  for  Tg/T^  *  10~^.  This  is  due  to  the  stabilizing  influence  of  the 
nonresonant  ion  contribution  (Eq.  (25)).  Also,  there  is  a  small  real 
frequency  associated  with  this  mode,  <  2  x  10  \  which  is  due  to  the 

density  gradient. 


The  weak  shear  case,  V  /0  L  =  -  0.01,  is  chosen  such  that 

o  e 


A  =  (k  /k  )  (k  /k  +  V  /S1L)  <  0  but  A  -  (k  /k  )(k  /k  +  V  /£2  L)  >  -  m  /ra. . 
zy  zy  oe  z  y  z  y  o  e  e  i 

We  anticipate  a  nonresonant  instability  in  the  "cold"  electron  regime  (Tg  « 
T^)  and  a  resonant  instability  in  the  "hot"  electron  regime  (Te  »  T^:).  This 
is  clearly  illlustrated  in  Fig.  3.  In  the  limit  Te  «  T^,  a  strong 

nonresonant  instability  occurs  which  asymptotes  to  a  growth  rate 

_3 

of  y/w  =  4.6  x  10  in  the  limit  T  +0.  Also,  there  is  a  small  real 

frequency  associated  with  this  mode  due  to  the  density  gradient  which,  for  the 

most  part,  has  c u  /a  <  1.0  x  10  ^  for  T_  <  T<.  As  T  /T.  1,  the  ions  become 

r  r  pe  e  i  e  i 

important  and  stabilize  the  mode  because  of  ion  Landau  damping, 


i.  e.  ,  a>  ~ 

kv . .  On  the 

i 

other 

hand. 

for  T  »  T,  the 
e  i 

resonant 

instability 

described 

in  Sec.  II.  B. 

2  is 

shown. 

Associated  with 

this  mode 

is  a  large 

real  frequency  (Eq.  (31)).  The  growth  rate  is  comparable  to,  but  smaller, 

than  the  nonresonant  "cold"  electron  mode.  Again,  as  T  /T.  *  1,  the 

e  i 

assumptions  used  to  arrive  at  Eq.  (32)  break  down  and  the  mode  is  stabilized 
by  ion  Landau  damping. 
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In  Fig.  4  we  plot  y  /u  vs.  k  r  for  V  /ft  L  ■  -  0. 01,  u  /ft  =  0.5, 
°  r  m  pe  y  Le  o  e  pe  e 


T  /T.  *  10  3  and  103 ,  and  r,  /L  -0  and  1.0  x  10  \ 
e  i  Le  n 


Here,  y  denotes  the 
m 


growth  rate  maximized  with  respect  to  k^/k^;  typically,  ~  0.005  for  the 

_3 

curves  shown.  We  first  discuss  the  "cold"  electron  case  (T  / T.  =  10  ). 

e  l 

Instability  extends  over  a  very  broad  range  in  k^  space  from 
-3 

k  r.  *  10  to  k  r,  =  1.  The  growth  rate  maximizes  for  k  r  *  0.1. 
y  Le  y  Le  °  y  Le 

As  k^r^  -*•  0,  the  instability  is  stabilized  because  of  ion  Landau  damping. 

The  additional  damping  due  to  a  density  gradient  is  evident  by  contrasting 

-4 

the  rT  /L  =0  and  r,  /L  =  1-0  x  10  curves  in  the  long  wavelength 

Le  n  Le  n 

limit  (k  rY  «  1).  The  mode  is  stabilized  in  the  short  wavelength 

y  Le 

limit  (k  r„  >  1)  because  of  finite  electron  Larmor  radius  effects, 
y  Le 

The  situation  for  the  resonant  instability  in  the  "hot"  electron 

3 

regime  (T  / T,  =  10  )  is  somewhat  different  In  this  case  the  waves  are  much 
e  i 

more  localized  in  k..  space ,  with  growth  mostly  occurring  in  the  narrow 

region  k  rT  *  0.1  -  1.0.  The  growth  rate  is  sharply  reduced  as  k  r  becomes 
y  Le  y  Le 

small  (i.e. ,  k  r.  <  0.1)  because  the  resonant  electron  contribution, 

7  ue 

proportional  to  exp(-u)2/k^  v^),  becomes  negligible  in  this  limit.  Again, 

for  k  r  >1,  the  mode  is  stabilized  because  of  finite  electron  Larmor  radius 
y  Le 

effects. 


In  Fig.  5  we  plot  y/u  vs.  v  / in  for  k  r  *  0.1,  k  /k  *  0.005, 

°  pe  e  pe  y  Le  z  y 

rY  /L  =  1.0  x  10"\  V  /ft  L  =  0.01,  T  /T.  =  10~3and  103,  and  u>  /ft  =0.5  and 
Len  oe  ei  pee 

2.0.  The  important  features  of  this  figure  are  the  following.  First, 
electron  collisions  eventually  stabilize  the  instability  regardless  of 
Te/T^.  Second,  electron  collisions  are  much  more  effective  at  stabilizing  the 
mode  when  “pe/^e  Is  large,  which  is  consistent  with  the  analytical  analysis 
presented  in  Sec.  II.  B-  1.  Finally,  there  appears  to  be  a  weak  collisional 

3 

instability  for  the  "hot"  electron  case  (T  /T.  =  10  )  when  o>  /ft  >1. 

'  e  1  pe  e 


Fig.  4  -  Plot  of  7m /wpe  vs,  kyr^  for  V0/neL  *  -  0.01,  =  0.5,  Te/Tj  =  10  3, 

103,  and  r^e/Ln  =  0, 1  X  10~4.  Note  narrowing  of  region  of  unstable  waves  for  Te/T; 
=  103. 


I 
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Fig.  5  —  Graph  of  vs.  J>e/cope  for  kyr^  =  0.1,  k2/ky  =  0.005,  r^/L,,  =  10~4, 
V0/J2eL  *  —  0.01,  Te/Tj  =  10“ 3 ,  103  ojpe/J2e  =  0.5,  2.  Note  reduction  in  growth 
rate  for  increasing  collision  frequency . 


III.  Summary  and  Discussion 


We  have  presented  a  linear,  electrostatic,  kinetic  theory  of  velocity 
sheared  electron  streams  flowing  parallel  to  a  magnetic  field.  For  simplicity 
wc  have  treated  the  case  where  the  electron  velocity  varies  transverse  to  the 
direction  of  the  flow.  Since  a  Vlasov  analysis  has  been  used,  finite  Larmor 
radius  effects  and  wave-particle  resonances  have  been  properly  treated.  These 
are  important  since  (1)  we  find  appreciable  wave  growth  for  kr  <  1,  (2) 

LG  ^ 

electron-wave  resonances  can  be  a  destabilizing  influence,  and  (3)  ion-wave 
resonances  are  a  stabilizing  influence.  In  addition,  we  have  included  an 
electron  density  gradient  transverse  to  the  flow  direction,  background  ion 
response,  and  electron  collisional  effects.  Moreover,  through  numerical 
solution  of  the  fundamental  dispersion  relation,  our  theory  is  valid  for 
arbitrary  values  of  Te/T^.  Our  principal  results  can  be  summarized  as 
follows: 

1.  For  "cold"  (Te  <<  T^)  velocity  sheared  electron  streams,  we  find  a 
nonresonant  instability  when 


3V 


if  ><0 


(33) 


We  note  that  a  density  gradient  is  stabilizing  and  provides  a  long 
wavelength  cutoff  such  that  the  unstable  waves  are  preferentially 
excited  locally  (k2L2  >  1). 


2.  For  "hot"  (Te  »  T^)  velocity  sheared  electron  streams,  our  results 
again  indicate  a  nonresonant  instability  when 
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(34) 


k 

k 


z 

y 


3V 


< 


3x 


m 

e 


m 


i 


As  in  Che  cold  electron  case,  the  density  gradient  exerts  a 
stabilizing  influence.  However,  we  have  also  found  a  resonant 
(electron-wave)  instability  that  has  the  sane  instability  condition  as 
the  cold  electron  case  (Eq.  (33)). 


3.  For  Te  »  ,  except  for  strong  shears,  (Eq.  (34)),  the  modes  are 

stabilized  due  to  ion  Landau  damping.  Ion  Landau  damping  also 
stabilizes  the  long  wavelength  modes  (i.e.  ,  k  r^e  0)  for 

T  /T.  *  0  . 
e  i 

4.  Electron  collisions  exert  a  stabilizing  influence  independent  of 

Tg/T^.  Collisions  are  more  effective  for  stabilization 

when  w  /2  >1. 

pe  e 

In  the  preceeding  theoretical  development  we  have  assumed  that  the 
unstable  modes  have  high  frequencies  y  >  such  that  the  background  ions 
cannot  execute  a  gyro  oscillation  on  the  instability  time  scale.  Thus,  the 
ions  are  assumed  unmagnetized.  For  y  <  the  ions  must  be  considered 
magnetized  and  the  previous  analysis  is  invalid.  In  this  case,  other 
instabilities,  e.g. ,  cyclotron  instabilities  [Kindel  and  Kennel,  1972],  are 
possible.  We  have  shown  a  posteriori  that  the  condition  y  >  0^  has  been  met 
for  the  parameter  regimes  studied. 

Nonlinear  numerical  studies  of  velocity  sheared  ions  flowing  parallel  to 
a  magnetic  field  have  recently  been  performed  (Tajima  and  Leboeuf ,  1980]. 

These  simulations  indicate  classical  vortex  formation  and  the  production  of 
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anomalous  viscosity  resulting  in  a  quasilinear  flattening  of  the  initially 
sheared  velocity  profile.  We  anticipate  similar  behavior  for  sheared  electron 
streams  and  defer  a  detailed  study  of  the  nonlinear  evolution  of  electron 
velocity  shear  instabilities  in  space  plasmas  to  a  later  report. 

For  application  of  this  theory  to  the  auroral  ionosphere,  we  consider  the 
interface,  usually  seen  near  discrete  arcs,  between  downward  flowing  high 
velocity  hot  electrons  and  upward  streaming  low  velocity  cold  return  current 
electrons.  Many  experimental  studies  of  auroral  arc  structure  and  dynamics 
have  been  made,  particularly  rocket  and  satellite  measurements.  Kelley  and 
Carlson  [1977]  have  detected  near  the  edge  of  an  F-region  auroral  arc, 
electrostatic  waves  with  spatial  scale  sizes  less  than  the  measured  ambient 
velocity  shear  scale  size.  These  results  are  not  inconsistent  with  our 
findings  of  a  local  (k2  L2  >>  1)  electrostatic  instability  due  to  transverse 
sheared  electron  streams  flowing  parallel  to  the  geomagnetic  field.  In 
addition,  several  other  investigators  [Whalen  and  McDiarraid.  1972;  Bryant  et 
al.  ,  1973;  Arno Icy  et  al.  ,  1974]  have  found  that  high  energy  field-aligned 

electrons  are  generally  localized  near  the  edges  of  auroral  arcs  in  regions 
separating  different  plasmas.  These  electrons  have  also  been  found  to  occur 
in  bursts.  Moreover,  electron  energy  spectra  near  the  edges  of  auroral  arcs 
[Bryant  et  al.  ,  1973  Bryant ,  1981  Carlson  and  Kelley,  1977]  often  show 

structure,  i.e. ,  cold  low  energy  electron  are  intermixed  with  hot  high  energy 
electrons  and  vice  versa.  These  features,  observed  near  the  edges  of  auroral 
arcs,  may  be  explained  by  the  linear  and  nonlinear  evolution  of  the  parallel 
electron  velocity  shear  instability.  This  can  be  seen  by  the  following 
scenario.  The  parallel  electron  velocity  shear  instability  can  evolve,  in  the 
nonlinear  regime,  into  vortices  which  can  be  described  as  a  turbulent  boundary 
layer  near  the  edges  of  the  electron  flow.  The  vortices  by  their  very 
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nature,  will  tend  to  mix  and  transport  hot  low  density  electrons  with  cold 
high  density  thermal  electrons.  These  vortical  structures  will  convert 
kinetic  streaming  energy  into  vortical  rotational  energy  and  thereby  act  as  a 
block  near  the  edges  of  the  electron  flow.  This  blockage  can  be  described  as 
a  "resistivity"  and  lead  to  the  formation  of  localized  electric  fields  and 
subsequent  acceleration.  As  the  vortices  grow  and  expand,  the  initial 
velocity  profile  will  flatten  in  a  quasilinear  fashion  due  to  anomalous 
viscosity  effects  [Tajima  and  Leboeuf,  1980;  iliura  and  Sato,  1978].  As  a 
results  this  "resistivity"  will  be  sporadic  and  bursty  in  character  as  the 
velocity  sheared  profile  alternately  steepens  and  flattens.  These  hypotheses 
will  be  tested  in  detail  in  future  work. 

Since,  for  conditions  typical  of  discrete  F-region  auroral  arcs,  i.e. 

-2  -2 

L  =  1-10  km,  k  /k  =10  ,  V  /ft  L  =  10  ,  we  find  unstable  wavelengths 
ranging  ranging  from  1=1-  100  m,  these  irregularities  could  be  observed 
using  radar  backscatter  methods  [Hanuise  et  al.  ,  1981].  We  note  that  since 
collisions  are  stabilizing,  this  electron  velocity  shear  instability  will  be 
preferentially  excited  at  high  altitudes.  In  addition,  throughout  our 
theoretical  treatment  of  this  instability,  we  have  ignored  beam  instabilities 
[Dungey  and  S  t rangeway  1976;  Papadopoulous  et  al.  ,  1974]  and  their 
interactions  with  velocity  gradient  driven  modes.  Finally,  these  results  may 
be  applicable  to  the  more  diffuse  interface  between  region  1  and  2  current, 
systems  in  the  auroral  ionosphere  [Iijima  and  Potemra,  1976]  and  to  electron 
current  return  current  regions  in  solar  flares  [Knight  and  Sturrock,  1977]. 
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01CY  ATTN  MARSHALL  CROSS 

HSS,  INC. 

2  ALFRED  CIRCLE 
3EDF0RD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
01CY  ATTN  J.M.  .AH IN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
50C  WASHINGTON  AVENUE 
NUTLEY,  SJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TCRREYAKA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 
01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAURAL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01  a  ATTN  THOMAS  POTEMRA 
•01CY  ATTN  JOHN  DASSOULAS 

Kama:;  sciences  corp 

P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 


01CY  ATTN  T.  MEAGHER 
KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


OICY 

ATTN  DASIAC 

OICY 

ATTN  TIM  STEPHANS 

OICY 

ATTN  WARREN  S.  KNAPP 

OICY 

ATTN  WILLIAM  MCNAMARA 

OICY 

ATTN  8.  GAMB1LL 

LINXABI 

T  CORP 

10453  ROSELLE 
SAN  DIEGO,  CA  92121 
01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO,,  INC 
P.O.  BOX  504 
SUNTIYVALE,  CA  94088 
0  ICY  ATTN  DEPT  60-12 
01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

C1CY  ATTN  MARTIN  WALT  DEPT  32-12 
01CY  ATTN  W.L,  LMHOF  DEPT  52-12 
01CY  ATTN  RICKARD  G.  JOHNSON  DEPT  5 
01CY  ATTN  J.3.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  S  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  4  ENGR.  CTR. 

4800  3 RADFORD  DRIVE 
HUNTSVILLE,  AL  358C7 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  SOX  5837 
ORLANDO,  FL  32805 
01CY  ATTN  R.  HEFFNER 

m.i.i.  li:;col:;  labohatory 

P.O.  30X  73 
LEXINGTON,  MA  02173 

0  ICY  ATTN  DAVID  M.  TOWLE 
0  ICY  ATTN  ?.  WALDRON 
01CY  ATTN  L.  LOUGKLIN 
OICY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  SOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 
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OiCY  ATTN  TECHNICAL  LIBRARY  SERVICES 
MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 

Santa  Barbara,  ca  93101 


0  ICY 

ATTN 

P.  FISCHER 

OICY 

ATTN 

U.F.  CREVIER 

OICY 

ATTN 

STEVEN  L.  GUTSCHE 

OICY 

ATTN 

D.  SAPPENFIELD 

OICY 

ATTN 

R.  30GUSCH 

OICY 

ATTN 

R.  HENDRICK 

OICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

DAVE  SOWLE 

OICY 

ATTN 

F.  FAJEN 

C1CY 

ATTN 

M.  SCHEI3E 

OICY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
?.0.  BOX  308 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORCANSTSRN 
0  ICY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WE5TGATE  RESEARCH  PARK 
1830  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
0 ICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1 A 56  CLOVERFIELD  3LVD. 

SANTA  MONICA,  CA  90404 
OICY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
218  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NC  CLASS  TO  THIS  .ADDRESS) 

•3 ICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHITOMETRICS,  INC. 

—  I  MARRETT  ROAD 
LEXINGTCN,  MA  02173 

OICY  ATTN  IRVING  L.  KCFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.0.  30X  3027 
BELLEVUE ,  WA  98009 
-02CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  30X  10367 
Oakland,  ca  946 10 
ATTN  A.  THOMSON 

310  ASSOCIATES 
P.O.  30X  9693 
MARINA  DEL  FEY,  CA  90291 
01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  8.  WRIGHT,  JR. 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
Q1CY  ATTN  C,  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
0  ICY  ATTN  CULLEN  CRAIN 
OtCY  ATTN  ED  8EDR0ZIAN 


RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
£0  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  AITS  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 
P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
0  ICY  ATTN  DANIEL  A.  HAMLIN 
0  ICY  ATTN  E.  FRISMAN 
0 ICY  ATTN  E.A.  STRAKER 
0 ICY  ATTN  CURTIS  A.  SMITH 
C1CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  12102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
.333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NEILSON 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

L.L.  C05B 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  G.  CKESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYS 

OICY 

ATTN 

M.  BARON 

OICY 

ATTN 

RAY  L.  LEAD  A3  RAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

J.  PETERSON 

C1CY 

ATTN 

R.  HAKE,  JR. 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
3EDF0RD,  MA  01730 
01CY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
SEDFCRD,  MA  01730 

01CY  ATTN  W.P.  3CQUIST 

TRW  DEFENSE  a  SPACE  SYS  GRCUF 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
01CY  AITN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VI SI DYNE 

SOUTH  3EDF0RD  STREET 
BURLINGTON,  MASS  01803 
3 ICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
0  ICY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  2037S 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
OR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  U.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  V.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CT  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92132 

MX.  R.  ROSE  -  CODE  S321 


NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  KISBET 
DR.  P.R.  ROHRBAUCH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08S40 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HELMERL 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
M.  PONGRATZ 
D.  SIMONS 
C.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 
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UNIVERSITY  OF  CALIFORNIA, 
LOS  ANCELES 
405  HIUGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIOSD1 
DR.  E.  OVERMAN 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  S CHUNK 
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